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Abstract 

We develop new more powerful techniques, based on an almost closed form for the 
lattice worldsheet propagator, for analyzing planar open string worldsheets defined on 
a lightcone lattice. We show that results obtained in earlier work are easily reproduced 
with far more precision. In particular, consistency checks which required numerical 
analysis in the earlier work can now be confirmed exactly. 
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1 Introduction 



The lightcone worldsheet [1-3] lattice was proposed long ago [4] as a method to digitize the 
summation of planar open string multiloop diagrams. Because the open string spectrum 
includes a massless spin 1 particle, this sum of diagrams should have the infrared behavior 
of large N [5] gauge theory. If the worldsheet lattice can reliably reproduce string theory 
diagrams, its a' — ?> limit should just as reliably reproduce gauge theory [6]. With this 
possibility in mind, we have recently embarked on a program [7] to critically evaluate the 
accuracy of this lattice in reproducing the continuum perturbative diagrams. In particular, 
it is very important that lattice artifacts be shown to be benign: they should either vanish 
in the continuum limit or be absorbed in redefinitions of parameters in the theory. This is a 
necessary prerequisite to applying these lattice methods to nonperturbative calculations of 
QCD. 

In [7] we studied the one loop self-energy diagram for the bosonic closed string in enough 
detail to see that the lattice accurately reproduced the ultraviolet behavior of the diagram. 
This is as much as we should expect, since the open bosonic string tachyon should and does 
ruin the infrared behavior of the diagram'^. The analysis in [7] employed what might be called 
a string field theory approach (see, for example [8]): the diagram was built up from open and 
closed string propagators. While this approach was manageable at one loop, it quickly gets 
unwieldy for multiloop diagrams. Even for the one loop open string self-energy diagram, 
attaining enough accTiracy to make definitive conclusions proved to be problematic. To 
improve on this situation we develop, in this article, a more powerful "worldsheet" approach 
based on the techniques of worldsheet quantum field theory defined on the lightcone lattice. 
The key to this approach is an almost closed form expression for the worldsheet propagator 
on the lattice (see Eq.(18)). 

The perturbative string field theory approach of [7] keeps manifest the contribution of 
all the intermediate string mass eigenstates contributing to the diagram. While the ensuing 
formulas for the self-energy shifts were exact at finite lattice spacing, we had to resort to 
numerical analysis to analyze the continuum limit. The extrapolation of our numerical results 
to the continuum was sufficiently accurate to make rather convincing consistency checks, such 
as a vanishing graviton self-energy. These checks were nonetheless subject to numerical error. 
In contrast, the methods of the present paper are powerful enough to analyze the continuum 
limit exactly in the ultraviolet and to confirm rigorously such consistency requirements. 

The Giles-Thorn (GT) discretization of the worldsheet [4] begins with a representation 
of the free closed or open string propagator as a lightcone worldsheet path integral defined 
on a lattice. The lattice replaces the transverse coordinates of the string x^a.r), living on 
a rectangular P"*" x T domain, with discretely labeled coordinates a;;^ = x{kaTQ, ja), living 
on an M X TV grid with spacing a, where = MoTq and T = a{N -\- 1). The free string 

"^It is logically possible that summing the bosonic string diagrams with an infrared cutoff stabilizes the 
vacuum in a way to produce QCD physics, but it is more likely that another string model, such as the 
tachyon-free Neveu-Schwarz sector of the superstring, is required to truly reproduce QCD. 
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propagator is then simply a Gaussian integral 




kj 

where the MN x MN matrix A is the lattice worldshect propagator that will be the central 
focus of this article. Then up to an overall normalization factor Dq — det"^^"^-* A~^, where 
D is the spacetime dimension [D = 26 for the bosonic string). 

On this lattice the sum of all open string multiloop planar diagrams can be obtained by 
summing over all patterns of missing spatial bonds. Formally, this is achieved by introducing 
Ising-like variables Sl — O,! and taking the worldsheet action to be 

^Planar = y E K^^'^' " ^^')' + "^^'(^m " ^/)'] 

+ {D - 2)B ^(1 -Si)~Y. [^.'(1 - S^'^') + "^.'^'(1 - ^9 (2) 

kj ij 

^ '^x^-[A-' + V{S)]x + A{{S}). (3) 

The terms in A{{S}) insert the coupling constant g in the appropriate way and allow for an 
open string self-energy counterterm B. Then we have 

V = voY,^^^~^\i +y^y~^^^^^^ ■ (4) 

{S} 

When K is a sparse matrix, i.e. when there are a relatively small number of missing bonds 
(Ylikji^ " ^kj) ^ M which can be arranged by taking B ^ 1), this will be a particularly 
efficient way to evaluate the terms of perturbation theory. Holding B sufficiently large serves 
as a physical and convenient infrared regulator in our studies of the properties of the planar 
diagrams. 

The paper is organized as follows. In Section 2 we construct the worldsheet propagator 
on the GT worldsheet lattice. It is remarkable that, in spite of the discretization of time 
{ix~^), the result is explicit and not much more complicated than the well known continuum 
worldsheet propagator. In section 3 we apply this expression to the calculation of the tachyon 
one loop closed string self energy. Our results, being exact, can be carried out to arbitrary 
precision, agreeing with the numerical results of [7] to the precision achieved in that reference 
(only 3 significant figures for some of the subleading contributions). The self-energy of the 
closed string graviton and selected higher mass states is similarly analyzed in Section 4. We 
conclude with discussion in Section 5 of the significance of our results and their promise for 
analyzing the open string self-energy as well as higher loop diagrams. 
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2 Lattice Worldsheet Propagators 



We develop the tools of quantum field theory for the worldsheet lattice. Of central interest 
are the worldsheet correlators of the coordinates on the M x N lattice corresponding to the 
free closed or open string, 



_ , , _ , _J±^-^i-^k ^ 



J Vx e 



(5) 



s 



Because the expectations are taken with Gaussian weight, the two point correlator in a 
single dimension captures all of the relevant information in arbitrary multi-point correlators 
in any number of dimensions. For the bosonic string we should of course take 26 space-time 
dimensions or 24 transverse dimensions. 

A straightforward evaluation is to use closure to write the numerator as the product of 
three string field propagators (see Appendix B): one from time —{N — j) to j, one from 
time j to /, and the last from time / to +(A^ + /). We can resolve x], x{ into normal 
modes g^, respectively. Then because each normal mode path integral is independent, 
{oinln) — ^mniQliQm) euds up with a Simple two variable Gaussian integral to do. 



/ 



{qLQu) = - j^2_B2 ^rnn- (6) 



Here A and B are read off from the formulas of Appendix B. For simplicity we set the g's 
at the initial and final times to zero. 
Then for non-zero modes they are: 

A ^ To sinh A [coth A^A + coth(Z - j)X\ , B = — -^^ — , (7) 

sinh(/ — jjA 

where A is A^ = 2 sinh^^ sin(m7r/2M) or A^ = 2 sinh~^ sin(m7r/M) for the open or closed 
string respectively. The non-zero mode contribution has a well defined N ^ oo limit: 

T sinh. A 

A ^ TosinhA[l + coth(/-j)A], B = , 

smh(Z — j)A 

-B If 1 A p-\^-i\'^ 



^2 _ ^2 sinh A V (cosh(Z - j) A + sinh(Z - j) A) ) 2Tq sinh A 

For the zero modes 

A - T^i±^ B ^° 
^0 — J-o^rrTi ^) -^0 



-Bo N (9) 

OT^ ATI '> \ ■> 



Al-B^ 2To{l + {l-j)/2N) 2To ATo 
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where we have taken A'^ large on the right side. To properly isolate the diverging term, we 
must remember that 2N is not the total time length of the string propagator. Rather the 
total length is 2Nt = 2N + l — Nx is the quantity that should be regarded as independent 
of j, I. In other words we should write 

Jj. i-i . ,.(^\ Nr 1-3 . I 



The zero mode contribution grows linearly with Nt- So it will be important that the zero 
mode be suppressed in the physical quantities that require the input of worldsheet propaga- 
tors. It is helpful to appreciate though that the divergent term is a constant independent of 
J, I. To define the inverse lattice Laplacian, it is consistent to drop it in effect modifying the 
boundary conditions on the Green function. 

In these derivations we have assumed I > j. For I < j the roles of the two indices are 
switched. When N ^ oo, we simply replace I — j ^ \^ ~ j\ the formulas: 

^-\l-j\Xm 

^^-^"^ ^ '^- 2TosinhA^ ' 

Prom their physical interpretation these are inverses of the lattice Laplacian^ 

(- A + 4 sinh^ A/2)/,- = 2/,- - - + 4fj sinh^ A/2 . (12) 

It is remarkable that this can be checked directly: 

r e-a-i)A(2-2coshA) l>j 
2e-l'-^l^ - e-l'+^-^'l^ - e-l'-^-^l^ - I e-(^-^)-^ (2 - 2 cosh A) l<j 

[2- e-^ - e-^ = (2 -2coshA) +2sinhA Z = j 

= -4e-l'-^lsinh2- + 25/jSinhA, (13) 

- A + 4 sinh^ - - — — = Si. , (14) 
2 J 2sinhA ' ^ ' 

which shows that (g^^'n) is the inverse of the lattice Laplacian on the nonzero modes. The 
proof for zero modes is even simpler 

-A\l-j\ = 2\l-j\-\l + l-j\-\l-l-j\ = -2Sij, (15) 

which confirms the same property for the zero modes. 



^Similar results for the inverse of the discrete Laplacian (discrete Green function) have actually appeared 
some time ago [9]. 
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Finally we return to the correlators on the spatial lattice by expanding in normal modes. 
The mode functions differ for the various types of string. For the Neumann open string 



M-l 



AO ml i i\ , , 2To v-^ , , ,, m(/i — l/2)7r mf/c — l/2)7r 

KiM = ^o«4)-^(gogo) + ^EtoL)cos M ''''^ M 

m=l 

Nt-\1- 3\ 1 e-l'-^l^™ mih - l/2)7r mik - l/2)7r 

m=l 



For the Dirichlet open string 



2Tn ^^^—^ , ■ , , . m/iTT . m/cvr 



M M M 

m=l 

M ^ smhA^ M M 

m=l '"• 

^mMi = ^o^^Ml/L) = ^— A^.,, = o, MM. (17) 



For the closed string 



^kki = WA) = + I Y.^A^mA'M-J exp '"^^^^ ^^^^ 

m=l 

= 2M +2M>^d^^"P M ■ 

m=l "* 

3 Closed String Self-Energy: Tachyon 

For the rest of the paper, we will apply our new approach in order to assess its calculational 
efficiency. In particular, we will use it to obtain 1-loop self-energy corrections to low-lying 
states of the closed string, for which we have a measure of comparison from our previous 
treatment [7]. In this section, we will particularly focus on the tachyon ground state. 



3.1 A Single Missing Link 

The matrix V has indices that are lattice locations, i.e. they are specified by two integers 
Vkj;mi- For a single missing link, at time j and linking spatial site k to site A; -|- 1, the term 
{To/2){xl._^i — is missing from 5". That means that 

ml;m'l' 
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from which we see 

yml;m'l' — —3lj3l'j{3m,k+l3m\k+l + 3m,kSm',k — Sm,k+l^m'k — Sm',k+l3mk)- (20) 

This matrix has entries only in rows and columns with labels kj and k + l,j, in other words 
a 2 X 2 submatrix. However the product matrix VA has nonzero entries only in rows with 
labels kj and k+l,j, but in general any column entry in these two rows can be nonzero: there 
are 2MN entries! But in calculating the determinant of / + V^A by expanding in minors, 
one quickly sees that it is only the 2x2 subblock of / + VA that contributes. Similarly if 
there are several missing links, the only part of / + VA that contributes to the determinant 
is a correspondingly sized subblock. 

Let us work out VA and the determinant for a single missing link. 



{VA)ml^qp — —Sij{Sm,k+lA(^k^i^j^qp + 5m,kAkj,qp — Sm,k+lAkj,qp — SmkA(^k^i)j^qp) 



— ~^lj{{^m,k+l — ^mk){^{k+l)j,qp 

Then the desired determinant is 



))■ 



det{I + VA) 

Prom (18) 

\k+i)j,kj 



det 
1 - A 



1 + ^{k+i)j,kj - ^kj,kj 



^kj,qp 



A(fc+l)j,(fc+l)j - Akj,(_k+l)j 



(21) 



^{k+l)j,kj 



+ A 



kj,kj 



^{k+l)j,{k+l)j 



kj,{k+l)j 



(k+l)j,(k+l)j 



M-1 



i + Akr 



kj,{k+l)j 



i + A 



{k+l)j,kj 



,■- A 



kj,kj ■ 



(22) 



^kj,kj 



—y 



2M ^ sinh 

m=l 



2mi7r 
exp — — 1 



M 



^ M-l 

--T 



sin^(m7r/M) 



M-l 



sinh 



2M 



M ^ 



sin(m7r/M) 



(23) 



m — - ^J'l + sin2(m7r/M) 

Evidently the same result is obtained for the difference Afej^(fc+i)j — A(fc+i)j^(fc+i)j. For large 
M we can apply the Euler-Maclaurin series 

^ M-l 



M 



TO=1 



\M 



^ (2A;)!M2fc 
where are the Bernoulli numbers, to get 



5:7Sl7T^(/'"-"(l)-/'"-"(0)). 



(24) 



^{k+l)j,kj 



^kj,kj 



1 




sin nx 




[/" 




2 


Jo 


/ 1 

V 1 + sin nx 


1 


'1 


n 


2 


2 " 


6M2 90M4 



2n 



+ 



Bo .3 

M6 ) 



B, 



24M4 



So finally 



det(/ + VA) = i + ^ 



+ 



TT 



90M4 



(25) 



(26) 
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Figure 1: GT worldsheet lattice for the closed string self-energy. The dotted 
lines are identified. There are K — 1 missing links, chosen for concreteness 
between the spatial positions k — 1 and k — M. 



3.2 Single Slit with K - 1 Missing Links 

The case of one missing link describes a one loop diagram with the loop occupying two time 
steps. A single loop occupying K time steps has K — 1 consecutive missing hnks, as we 
depict in figure 1. Proceeding with the case of K — 1 missing links, again between spatial 
sites k and k+1, but this time for the time interval between instants J + 1 to J + — 1, it is 
evident that we will simply have to sum the right-hand side of (20) over j G + J+K — 1]. 
The sum will also carry over to (21), whose nontrivial subblock will now have size 2(K — 1). 
Due to the difference of delta functions on the latter relation, clearly the matrix rows with 
m ~ k + 1 will have the opposite values of the rows with m — k. So sorting our rows such 
that 

(m, I) = {{k, J + 1), . . . (A;, J + X - 1), . . . , (A; + 1, J + 1), . . . , (A;, J + X - 1)} , (27) 
and similarly for the columns, we can write in [K — 1) x [K — 1) block form 

det(/ + l/A) = det(^+/ ,5^)=det(/^ ^ 

= de'(_A / + ^-s) =''''*'^ + '*"-^'' 

where we employed elementary row and column manipulations that leave the determinant 
invariant, together with the block matrix identity 

det ( ^ M = det(g) det(5) . (29) 
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K 


det(/iip(a;)) 


2 
3 
4 
5 
6 


_2- 64+ 

16 ?f + 3' + ( 64 7> + + ^3^2) X 

IOQ 1441792 1 4I056 1 5o72 '^512 , 7 7«q 83886S8 , 26'2l44 , 1^3840 1024 ^ _ 
8l7r4 1 277r3 ' 7r2 37r ' I ' SItt* I 277r3 ' 97r2 tt j 



Table 1: Asymptotic expansion up to 0{x), x = of the determinant (36), 

entering the tachyon energy shift (38), for a sht of length if — 1, K = 2, ... ,6. 
Evidently, the coefficients of the expansion can be calculated exactly. 



In formula (28) above, 

— ^{k+l)l,kp — ^kl,kp , Bip = A(fc+i);^(fc+i)p — Akl^(k+l)p , (30) 

and since these quantities depend on /,p only through |/ — p\, the value of J + 1 will be 
immaterial and we can set it to zero. Hence our final expression for the determinant will be 

det(/ + 1/A) -det(/i,p), /,p = 1, 2, . . . K - 1 , (31) 

where 

hip — ^Ip + A{k-\-l)l,kp — Akl,kp + Akl,{k+l)p — A[k+l)l,{k+l)p ■ (32) 

Notice in particular that we now have the determinant of a [K — l)-dimensional matrix, 
whose elements depend on differences of propagators, such that the zero modes in (16)-(18) 
always cancel out. Specializing to the case of the closed string propagator (18), it is easy to 
show that^ 

= - ^ E y/'^'^T/^L, f sin(m7r/M) + ^1 + sin^(m7r/M)) " , (33) 
yl + sin^(m7r/M) V / 

and applying again the Euler-Maclaurin formula (24), we obtain 

(-l + 3|/-pn7r3 l\ ,,,, 



where we have reduced the integral I\i-p\ to a simple finite sum in appendix C. Reexpressing 

TT 

6M2 



hip{x) = hip{Q) + X + O [x^) , ^ = , (35) 



we can separate the M-dependence of the determinant. 



dei{hi,) = detihi.m + detih,,ix))l^, + O (^^) . (36) 



^Evidently, hip = h{\l —p\), namely det{hip) is the determinant of symmetric Toeplitz matrix. It can be 
shown [10] that it further reduces to a product of two determinants of approximately half the size. 
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K 


^-^Cdosed 


-^Pcdosed actual 


2 
3 
4 
5 
6 


0.1044844648 - 1.31291/Af2 
0.027700432 - 0.9578/Af2 
0.010959556 - 0.7268/M2 
0.005388196 - 0.5811/M2 
0.003032942 - 0.4828/M2 


0.104484465146 - 1.31299/M2 
0.0277004334342 - 0.957933/M2 
0.0109595576932 - 0.727031/M2 
0.00538819758183 - 0.581471/M2 
0.00303294412639 - 0.483277/M2 



Table 2: Asymptotic expansion up to 0(1/M^), of the tachyon energy shift sum- 
mand (38) for K = 2, ... ,6, where K — 1 is the slit length. The LHS coefficients 
were determined by fitting the M dependence, as in [7], with an error estimate at 
the order of the last digit. The RHS coefficients have been calculated exactly with 
the methods of the present paper, and evaluated up to the desired precision. 



Finally, the summand of the energy shift will be given by 

^-2A{K-l)Bo g-24(ii--l)Bo 



dee\hi,{0)) [ M2 det{hi,{0)) J ^ \M^J ' ^ ' 

where hip is given exactly in (33) and asymptotically in (34), and similarly we have expressed 
the summand in its exact (37) and asymptotic (38) form. 

Let us now discuss how our current, worldsheet-based approach, compares to the string 
field theory-related approach we employed in [7], when it comes to computing the tachyon 
energy shift. Clearly, the formulas wc have derived here involve determinants of size roughly 
K, so that they are advantageous for analyzing the summand in the ultraviolet region K <^ 
M. Conversely, the approach [7] yields determinants of size M, more suitable for the infrared 
K ^ M regime. 

The fact that each approach is more suitable for one of the two domains, is also evident 
in our ability to derive asymptotic formulas there. Because in our earlier paper we had no 
such formula for analyzing the K <^ M behavior of the integrand, we had to use the exact 
formula for the summand, evaluate it for a range of M and and perform fits in both 
variables in order to find the respective dependence. In contrast, here we obtain explicitly 
the form of the asymptotic expansion in M, and we only need to fit for the dependence of 
the coefficients on K. 

Furthermore, it is possible to compute these coefficients exactly for specific iiT, or evaluate 
them at arbitrary precision. As an example, we present the exact values for the first two 
coefficients of det(/i;p(x)) for X = 2, . . . , 6 in table 1. In table 2, we also present the respective 

asymptotic expansion for the summand (38), and compare the values for the coefficients, on 
the one hand obtained with our current method, and on the other hand by fitting the M 
dependence along the lines of [7]^. 



^See also figure 5 in the latter reference. 
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Clearly, the two results agree within our margins of error, notice however that their 

difference increases with K. This is a result of the systematic error coming from not taking 
into account the 0{1/M^) term in the fits, whose relative size also increases with K. Apart 
from the asymptotic expansion, we also confirmed that the exact formulas for the summand 
of the tachyon energy shift in the two approaches, (37) here and (51) in [7], agree for a large 
set of M, K values. 



4 Closed String Self-Energy: Graviton 

To extract information, e.g. energy shifts, about excited closed string states, we will need 
to consider the propagator on a worldsheet which includes interactions. If we denote the 
matrix describing a particular configuration of missing links by V , as we did in (3), then the 
propagator in question will be given by 

= (A-i + \/)-^ = A(/ + \/A)-i = A-A(/ + yA)-VA = A- AVA. (39) 

The final form on the right is useful when V is sparse, because then the inverse matrix 
appearing in the second term can be evaluated as the inverse of the submatrix obtained by 
projecting onto the sparse subspace. 

Using index notation for the propagator, A^^^, and choosing the times q and j much 
earlier and much later than all of the times occupied by V respectively, we can also write 



A^.. - ^ exp |-jA„ + .A„, + -(n,^ - Wrfj ■ (40) 

m,m 

2M^sinh sinh A^' 

Vmm' = ^ exp I IXm " sX^' " ^("^^ " "^'^) } ^kl,rs ■ (42) 
kl,rs 

Wc have normalized Aj^„j, so that it is 5mm' for V = 0. Thus A^^, det~''"^(/ + V^) gives the 
probability amplitude that the mode m' at early times evolves to mode m at late times. For 
the graviton self-energy, the relevant process is modes m = 1, M — 1 at early times evolving 
to the same modes at late times. Thus this contribution to the graviton self energy is 

- (ai:;a;;,_,)(^_,) + Ar(M-i) a^m-di) + ^a) . (43) 

4.1 A Single Missing Link 

For starters, let's take V with a single missing link. Its nonvanishing 2x2 subblock is the 
matrix 
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Putting A = A(k+i)3,k3-^kj,kj = ^{k+i)j,kj- ^{k-^i)j,{k+i)j the matrix 7 + FA projected onto 
the subspace of V and its inverse times V are 



/ + 1/A 
V = (/ + 1/A)-V 
Then we easily compute 



l + A 
-A 

1 

1 + 2A 



-A 
1 + A 

1 + A A 
A 1 + A 



-1 1 
1 -1 



V 



1 + 2A 



(45) 



^^:i(Xm-x^,)+ni(m'-m)(2k+i)/M ^HT^rn/M) sm{7rm'/M) 

1 + 2A 

sin^(7rm/M) ~ _ ,^_2^im(2fc+i)/M sin^(W^) 

l^m(M-m) — 1 + 2A 



-4 



1 + 2A 



(46) 
(47) 



Then we find 



V 

(M-m)(M-m) 



1 + 4 
1 + 
-4 



sin^(7rm/M) 
2M(1 + 2A) sinh A, 
27rm 



M2 3M4 
sin^(7rm/M) 



2M(1 + 2A) sinh A, 
27rm 



-2-Kim{2k+l)/M 



(48) 



M2 



27rim(2fc+l)/M , q ( ^ \ 

) 



(49) 



The one missing hnk contribution to the self energy of the closed string state |m, M — m) ± 
\M - m, m) is up to 0(1/M^) 

^mm^(M-m)(M-m) =^ ^m(M-m)^(M-m)m 



dct'^(l + V/^) 
(1 + inm/M^ - 4m7r2 (1 - 6m + 2m27r)/3M4) 
(1/2 + 7r/6M2 + 7r3 /90M^)^^ 

47r(m - 1) 27r2(-65 + 130m - 30(1 ± l)m2 + 2n + 2Qm^n] 



-2^2 I 1 + 



M2 



15M4 



.(50) 



This formula includes the shift for tachyon (m = 0) and the graviton (m = 1). The latter 
receives no l/M^ correction, consistent with zero shift in the continuum limit. Note also 
that we have assumed in these formulas that the polarizations of the first and second entries 
of |m, m') are different, so they don't properly describe the dilaton self-energy shift. 

Furthermore, for the graviton (m = 1 and plus sign), we can also compare the above 
formula with the fits we obtained for the value of the graviton energy shift in [7] . Multiplying 
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(50) with the boundary counterterm exp(— 24So), and writing the result in the notation of 
the latter paper, we have 

The result above is in excellent agreement with the fit presented in figure 11 of [7]. 
4.2 Single Slit with K - 1 Missing Links 

Let us now generalize the discussion of the previous section, for the worldsheet configuration 
where a link is missing between the same spatial sites k and A; + 1, but for a time interval 
K — 1 links long. Using the same reasoning as for the tachyon in the same configuration, it 
is possible to show that the matrix V defined in (39) has the special structure 

Va;;,A;s = V(fc+i);,(fc+i)s = —V{k+l)l,ks = —^k+l,{k+l)s = , (52) 

where h is the same {K — l)-dimensional matrix appearing in (31). With the help of these 
relations, we do the k, r summation in (42) 

V^^, = _ ^ e'-^'""^^-' ^e^^'"'""')'' + e^K-He^+l) _ ^^{.m'+m'k-mk) _ ^^{m'-mk-m)^^ f^^l^ 
Is 

and if we take out an overall factor exp[7ri(m' — m){2k + 1)], this simplifies to 

V^^, = _4e-K— sin ^ sin ^ ^ ^i>^m-sx^, j^-i (53) 

Is 

It is evident that the above relation implies 

—m){M—m) = Vmm , (54) 

so that our final working formula for the graviton summand, also including the required 
boundary counterterm, will be 

-24(K-l)Bo . X 

xp- - __1 (A^A^ +A^ A^ 1 

^-nSraviton - , 12/j , a n ^^11^(M-1)(M-1) + ^1(M-1)^(M-1)1 J 



dee\l + VA) 

g-24{A--l)So 

det^^(/ + \^A) 

.-24(ft:-l)So 



1 



2 / ~ \ 2' 
Vii \ , / Vii 



2MsinhAiy V2MsinhAi 
12. r , T^A^ (l + 2t/ + 2t/2) , 



(55) 



det^\l + VA) 
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K 


^-^Graviton 


"■'^-Pcraviton actual 


2 
3 
4 
5 
6 


0.104484465145- lO.ig/M^ 
0.027700433434 - 3.85/Af^ 
0.010959557693+ 1.87/M^ 
0.0053881975+ 6.82/M4 
0.003032944127+ 11.28/M^ 


0.10448446514630- 10.1905/M4 
0.02770043343416 - 3.8499/Af^ 
0.01095955769317+ 1.8837/M'i 
0.00538819758183 + 6.8571/M4 
0.00303294412639+ 11.3355/M4 



Table 3: Asymptotic expansion up to ©(l/M^), of the graviton energy shift sum- 
mand for X = 2, . . . , 6, where K — 1 \s the sht length. The LHS coefficients have 
been determined by fitting M, as in [7], with an error estimate at the order of the 
last digit. The RHS coefficients have been calculated exactly with the methods of 
the present paper, and evaluated with two additional digits of precision. 



where 

• TT K-l / 

sm X? / 

11 = ■ M V ( .sin , . , 

M V M 



■ TT K-l / I s 2{l-s) 

sm X? / . TT I . o TT 



MWl + sm^^ I,. 



and hr^ is the inverse of the {K — l)-dimensional matrix with elements (32). The same 
interesting phenomenon that we encountered for the tachyon also appears here, namely we 
can reduce the size of the matrices entering the energy shift by a half. The asymptotic 
expansion of (55) in M readily follows from the respective expansion of his (34), and in 
particular it is easy to show that 

= V.'(o) - ^ (e "sHo)) ".'(0)) +o{±-^, (57) 

where /i(0) is the M-independent part of the matrix h. Because the overall factor in (56) 
starts as (9(1/M^), we do not need additional terms in order to obtain (55) at 0{1/M'^). 
In fact, if we only focus at 0(1/M^) for a moment, the term on the right-hand side of (55) 
simplifies to 

1 + 2[/ + 2[/2 ~ 1 + — V h-HQ) - 1 + idet(V(x))|^^o 

1 + 2ty + 2ty - 1 + ^2-. ''is (U) - 1 + ^2 det(/i,,(o)) ' ^^^^ 

where for the last equality we used the identity 

-det(/.) = det(/.)TrU-i- , (59) 



dx \ dx ^ 

and also the fact that in our case the derivative matrix has all entries equal to one. 

Comparing (58) and (55) with (38), we observe that we have rigorously proven two 

important facts: That the leading order of the asymptotic M expansion for the graviton is 
equal to the tachyon one, and that the subleading term is always 0{1/M*) for any K -C 
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M. Of course these properties were expected to hold on physical grounds, however in the 
approach of our previous paper, we could only obtain empirical indications about them from 
the fits. Finally, for sample slit lengths, we compare the coefficients of the aforementioned 
fits with the exact values obtained with our new method, and evaluated at higher precision, 
in table 3. 

5 Discussion and Conclusion 

In this paper, wc continued our investigation of lattice-regularized string theory in the light- 
cone gauge, by introducing a new approach for evaluating the corresponding path integral. 
Whereas in our earlier work [7] we built the path integral by integrating products of free 
string propagators over the interaction points, here we treated it as a quantum field theory on 
the worldsheet. Given that free string propagators are the two-point functions of string field 
theory, we could call the former approach string field theory-based, and the latter approach 
worldshcct-bascd. 

The key idea for treating string interactions in this framework, was to examine how the 
path integral is modified as we start removing finks from the free worldsheet (3) . An essential 
ingredient for describing this departure, is the worldsheet correlation function A of two target 
space coordinates (5). We consider as the main result of this paper, the determination of 
this quantity explicitly in Fourier mode space (11), and as a simple sum in coordinate space 
(16)-(18). 

We then moved on to assess the efficiency of the worldsheet approach, by calculating the 
one-loop self-energy corrections for the closed string tachyon and graviton, and performing 
a comparison with the results of [7]. The self-energy corrections involve determinants of size 
M and K for the string field theory and worldsheet approach respectively, and hence the 
first one is more convenient for analyzing the K ^ M (infrared) regime, whereas the second 
one for the K <^ M (ultraviolet) regime.^ 

Indeed, with our current approach we were able to not only find the structure of the 
asymptotic expansion in M of the self-energy summand, but also to calculate its coefficients 
exactly, for each value of K. In this manner, we were able to rigorously prove two important 
facts, for which we only had strong indications up to now. Namely, that the leading term in 
the expansion is the same for the tachyon and the graviton, and that the 0{1/M^) subleading 
term for the graviton is zero (i.e. the graviton is massless in the ultraviolet region). In 
contradistinction, analyzing the K <^ M regime in [7] had to rely on fits for both variables 
M and K, which introduced larger numerical errors and made conclusions less definitive. 

Apart from its calculational virtues, our new approach adopts the point of view, implicit 
in our representation of the planar sum as a sum over Ising spin variables (2), which is 
much closer to the treatment of more general lattice systems: Each string diagram is like 
a lattice state, and we build all states by gradually adding more and more 'excitations', 
namely missing links, to the 'vacuum', or free worldsheet. It would be very interesting to 

''We remind the reader that M is the spatial size of the lattice and K is the temporal length of the slit 
representing the loop. 
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explore this point of view further, as it seems to suggest that string diagrams of different 
loop order but same excitation number may be similar to each other. Indeed, generalizing 
the considerations of sections (3.2) and (4.2) for arbitrary positions of the K — 1 missing 
links (m, /) = {(fci, ji), . . . , {kx-i, Jk-i)}, yields again a {K — l)-dimensional determinant, 
this time with elements 

hml,m'l' — Smm'^mm' + ^{m+l)l,m'l' ~ ^ml,m'l' + ^ml,{m'+l)l' ~ ^{m+l)l,{m'+l)l' ■ (60) 

It may be more advantageous to organize the sum over diagrams not by loop order, as 
dictated by the conventional wisdom of string perturbation theory, but by the number of 
missing link 'excitations'. 

With this more efficient method now in place, a primary objective will be its application 
for the study of the one-loop self-energy corrections to the low-lying states of open string 
theory, which though more intricate, is of main interest because of its relation to large N 
gauge theory. Once wc have similarly established the compatibility of the lattice regular- 
ization with Lorentz invariance in this well, then the next natural step will be the 

numerical evaluation of the full path integral with the help of Monte Carlo methods. 

In this respect, it will be very interesting to examine whether efficiency can be further 
improved by performing the sums in (16)-(18) analytically, in order to obtain explicit ex- 
pressions for the worldsheet propagators in coordinate space as well. In the most probable 
scenario, that the summation of all bosonic string diagrams does not succeed in stabilizing 
the vacuum, we will of course be aiming to develop a similar treatment for the superstring 
as well. 

Acknowledgements: This research was supported in part by the Department of Energy 
under Grant No. DE-FG02-97ER-41029. 

A Normal Modes 

A string with P+ = MqTq is described at a fixed time by M coordinates Xi or yi, i = 1, • • • M. 
In this article we require several normal mode decompositions depending on the boundary 
conditions. 

Neumann Open String 



1 



5o + \/ T7 




m=l 



mTr{i - 1/2) 
M 



(61) 




(62) 
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2 '^—^ . rmrk 



Dirichlet Open String 

Vk = Vm^^d^ ~Jf A; = 1, • • • , M - 1, yM = Qdm , (63) 

m=l 

gDm = V M ^ < m < M, Qdm^Vm- (64) 
fc=i 

Closed String 

1 '^-^^ 2mk'rii 1 ■r^^ 2(^1/ — n/jkni 

- ^l^^rne^v^^, yl^ = -^^a;.exp . (65) 

' m=0 ^ k 

This goes to the normal mode expansion with trigonometric functions with the substitutions 
Am = A*M-m = i^cm " iqsm)/V2, with < m < M/2, Ao = go, and Am/2 = gcM/2 (if M is 
even) . Prom this dictionary the non-zero correlators are 

1 M 

{A„AM-m) = -{{qcrr^qcn.) + {qsrnqsn.)) ^ ^ 0, - , 

{AqAo) = {qoqo), {AM/2AM/2) ^ {qM/2qM/2) ■ (66) 

B String Propagators 

B.l Neumann Open String Propagator 

(iV + l,a;^|0,xT'''" = popen^^ ^ ;L)giiy„pe„ ^ (57) 



w — -^0 

^ open — 2 
M-1 



(go,/ - 90,0 



2 



+ E ^i-h a;; ((g^, + qlj) coth(Ar + i)a;;, - ] , (68) 

m=l ^ V / m/ 

XI = 2sinh-i(sin— ) . (69) 

Where the q„iS arc the normal mode coordinates for the x's. The right side is the result of 
doing the integrations over all the xl with i = 1, ■ ■ ■ , M and j = 1, - ■ ■ N. The propagator 
spans N + 1 time steps and this result corresponds to assigning half the potential energy 
To Y.Zi\4+i - 4)y^ to time j = and half to j = + 1. 

B.2 Dirichlet Open String Propagator 

The Dirichlet open string propagator over a time of K = N + 1 steps is evaluated to be 

{qf,N + l\q\0)'' = V''{N+l)e'^'' , (70) 
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where 



27ry 



m=l 



sinh A5 



(71) 



(72) 



.-1 1 



2sinh-^-==, A^ = A^ = 2sinh-^sin— , m = 1, • • • M - 1 . (73) 



We recall that the above expressions give the the result of integrating over all the variables 
ul^ for j = 1,---,N, with half the potential energy assigned to j — 0,N + 1, which is 
consistent with the closure requirement. 

B.3 Closed String Propagator 



inclosed — 
M-1 



To 
' 2 



f - 



N + 1 



ivi—i / 

+ E ^^^^ ( (^-.^ + ^-./) ^°*^(^ + l)^m - 2 
m=l ^ 



Qm,iQm,f 



smh{N + 1)A 



A^ = 2 sinh ^ (sin^) 



(74) 

,(75) 
(76) 



Where the ^^'s are the normal mode coordinates for the x's. When we divide the closed 
string normal modes into sine and cosine modes, we arbitrarily call the m > M/2 modes 
sine modes and the m < M/2 modes cosine modes. When M is even, the M/2 mode is not 
doubles. The right side is the result of doing the integrations over all the with i = 1, ■ ■ ■ , M 
and j = 1, - ■ • N. The propagator spans + 1 time steps and this result corresponds to 
assigning half the potential energy Tq X]i^i(^i+i ~^i )^/2 to time j — and half to j = A^+1. 
In sums like these it is understood that x-^j^j_^_-^ = x\. Whenever we concatenate at a time j 
propagators with different numbers of missing links, we will understand that we add terms 
To(Ax)^/4 in the exponent so that the potential assigned to time j is that of the system 
with the least number of missing links. For example, the concatenation of an open string 
propagator with a closed string propagator entails the addition of Tq^x^^^ — x{Y /A to the 
exponent. 



C A Useful Euler-Maclaurin Expansion 

As we saw in in section 3.2 for the closed string propagator, when many links are missing, 
the leading term in the Euler-Maclaurin expansion of the elements of the corresponding 
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n 


In 





1/2 ~ 0.5 


1 


-l/2 + 2/7r ~ 0.1366 


2 


-5/2 + 8/n ~ 0.04648 


3 


-25 /2 + 118/(37r) ~ 0.02019 


4 


-129/2 + 608/(37r) ~ 0.01080 


5 


-681/2 + 16046/(157r) ~ 0.006696 


6 


-3653/2 + 86072/(157r) ~ 0.004568 



Table 4: Values of integral /„ (77), forn = 0, 1, . . . , 6. 



determinant involves an integral of the form 



, sinnx / . r \-2n 

In = ax I sm TTX + V 1 + sm irxj [77) 

Jo V 1 + sin TTX ^ ' 

f\ sin f / TTx r ~2 
— / ax — . = I sm h W 1 + sm 

•^0 ^i + sin^i^V 2 V 

/■\ Sin^ / . Tlx r . r,TTx\' 

= / dx = -sm — + Wl + sm^ — 

io ^i + sin^f V 2 V 2; 

_ 2 



-2n 

TTX 



2n 



We can evaluate this with the help of the identity [10] 



2^ n n 
+ V^TT^) " = Yl ^r^rZ''' + V^TT^ J] ^lnrZ'''-^ , (78) 



r=0 r=l 

where 

n f n + r \ 2r 



A„. = -- ' 2^^ (79) 



so that the integral can be rewritten as 



- Z.^-20Fr(l + i) Z^^-0Fr(i + r)- ^''^ 

If desired, we can formally express these finite sums in terms of hypergeometric functions, 
for example 

'"-^^^TTT^ = n 2Fi(l - n, 1 + n; 2; -1) . (81) 
^ V7rr(l + r) 
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In any case, the sums can be readily evaluated for specific values of n, and for the reader's 

convenience we have tabulated the first few cases in table 4. 

Summarizing, the sums that are relevant for the computation of the closed string propa- 
gator when many missing links are present, have an Euler-Maclaurin expansion of the form 
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